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Abstract 

We investigate, both experimentally and theoretically, the bifurcation to alternans in heart 
tissue. Previously, this phenomenon has been modeled either as a smooth or as border-collision 
period-doubling bifurcation. Using a new experimental technique, we find a hybrid behavior: very 
close to the bifurcation point the dynamics are smooth-like, whereas further away they are border- 
collision-like. This behavior is captured by a new type of model, called an unfolded border-collision 
bifurcation. 

PACS numbers: 87.19.Hh, 05.45.-a, 87.10.+e 
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Many nonlinear systems display a bifurcation, where the system's response changes qual- 
itatively as an adjustable parameter - the bifurcation parameter - is varied [lj. Bifurcation 
theory provides useful tools to understand the behavior of the system in the vicinity of the 
bifurcation. The most widely investigated nonlinear systems are those described by smooth 
differential equations or maps. As a result, there exists a vast literature on the classes of 
bifurcations occurring in these systems. On the other hand, bifurcations occurring in piece- 
wise smooth systems, such as border-collision bifurcations, are still under investigation and 
are known to display much richer behaviors 2]. Border-collision bifurcations are believed 
to occur in a variety of systems, such as mechanical and electrical devices that involve sud- 
den switching of a component, and in economic systems where decisions are made when a 
variable crosses a threshold. 

The primary purpose of this Letter is to describe experimental observations and a re- 
sulting new model of the bifurcation to alternans (defined below) in paced bullfrog cardiac 
muscle. We find that the bifurcation mediating this transition displays both smooth and 
piecewise-smooth characteristics. Specifically, by investigating the system's sensitivity to 
perturbations, we find smooth-like behavior near the bifurcation point while, further away, 
the behavior is border-collision like. This apparently contradictory behavior is placed in 
the context of a new type of model, which we call an unfolded border-collision bifurcation. 
The bifurcation to alternans is a crucial "problem to study because there is evidence that it 
can initiate ventricular fibrillation DDE! , which often underlies sudden cardiac death, one 
of the leading causes of death in the United States jf|. Therefore, determining the type of 
bifurcation is a crucial step in reaching the ultimate goal of suppressing cardiac alternans 



7] 



Before describing our findings, we review briefly the behavior of paced cardiac muscle. An 
applied electrical stimulus induces an action potential, which is the characteristic time course 
of the transmembrane voltage. In experiments, the dynamical state of the tissue is often 
described by measuring the action potential duration (APD). The pacing interval, known 
as the basic cycle length (B), serves as the bifurcation parameter. Under various conditions 
of periodic electrical pacing, M stimuli can elicit N responses (M:N behavior) for different 
values of B, where the transition from one response pattern to another under changes in B 
is mediated by a bifurcation. For slow pacing, 1:1 behavior is usually observed, where each 
action potential is identical. As B shortens, a period-doubling bifurcation sometimes occurs 
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FIG. 1: Schematic bifurcation diagrams with discrete sampling. The sampled points (solid dots) 
are identical for the (a) smooth and (b) border-collision bifurcation, (c-f) Alternate pacing: the 
trend in V vs. B is illustrated for a (c) smooth and (d) border-collision bifurcation and the trend 
in r vs. 5 is illustrated for a (e) smooth and (f) border-collision bifurcation. 

8], giving rise to a long-short periodic pattern in APD known as a 2:2 rhythm or alternans 

Si- 
Beginning with the pioneering work of Nolasco and Dahlen jlol |. it is generally believed 
that the transition to alternans is due to a smooth supercritical period-doubling bifurcation. 
More recently, Sun et al. [llj introduced a model for atrioventricular nodal conduction 
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in which it has been shown [7|, [12j that the transition to alternans arises from a border- 
collision bifurcation. The bifurcation diagrams, where the steady-state values of APD are 
plotted as a function of B, are different for these two cases. For a smooth period-doubling 
bifurcation, the bifurcated curves become tangent at the bifurcation point, as shown in Fig. 
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1(a), whereas they are generally not tangent for a border-collision bifurcation, as shown in 
Fig. 1(b). Unfortunately, it is exceedingly difficult to measure the bifurcation diagram with 
sufficient resolution to distinguish between these behaviors experimentally because only a 
finite number of pacing intervals can be tested before tissue damage occurs. For example, 
Figs, ^a) and (b) show that the same set of data (solid dots) is consistent with both types 
of models. 

We apply a new technique involving alternate pacing to investigate the bifurcation 
type in cardiac tissue. This method is a modification of a technique introduced by Heldstab 
et al. for general dynamical systems. Specifically, this method relies on periodic per- 
turbations to the bifurcation parameter to investigate the dynamics of the system in the 1:1 
region. More recently, Karma and Shiferaw suggest that Heldstab's method ^| can be 
exploited in a clinical setting to study the likelihood of an unfavorable cardiac event. 

Our alternate pacing protocol is implemented in the following way. In the 1:1 regime, 
we perturb the nominal pacing interval Bq by a small value 8 so that, for the n th stimulus, 

B n = B + (-l) n 5. (1) 

As a result, APD alternates in a long-short pattern with the corresponding steady-state val- 
ues denoted by APD\ ong and APD^^. To measure the system's sensitivity to perturbations, 
we define a gain as 

~ 25 ' U 

In previous papers |l3t llq . we have shown that, for the two types of bifurcations, T 
depends on the bifurcation parameter B and perturbation amplitude 5 in qualitatively 
different ways. Specifically, we computed T under variations of both B and 5. For a 
smooth bifurcation, T decreases monotonically as Bq increases compared to i?bif [see Fig. 
mc)]. This is because the bifurcation parameter's sensitivity to perturbations decreases 
further away from the bifurcation point. Likewise, V decreases monotonically as 5 increases 
[see Fig. life)]. Specifically, T is proportional to 5~ 2 ' 3 for a range of 5, and it saturates 
for sufficiently small 5 [16]. On the other hand, for a border-collision bifurcation, T is a 
piecewise smooth function of B and 5. Overall, when the alternating response crosses the 
border (which occurs at B cr u and 5 cr u), T exhibits a decreasing trend as Bo increases [see 
Fig. CJd)] but an increasing trend as 5 increases [see Fig. dT)] [l^]. Although Figs. [T^c) 
and (d) are slightly different, it is difficult to distinguish their differences with discretely- 



sampled data. However, because of significant dissimilarity between Figs. |IJe) and (f), the 
distinction between the two bifurcation types is evident even with just a few data points 
and in the presence of experimental noise. 

We apply the alternate pacing protocol in 6 adult bullfrogs. In these experiments, the 
heart is excised from an adult bullfrog (Rana catesbeiana) of either sex. After the pacemaker 
cells are cut away, the top half of the ventricle is removed and placed in a chamber that is 



superfused with a recirculated physiological solution 



A bipolar extracellular electrode 



applies 2-ms-long rectangular current pulses to the epicardial surface of the tissue [3| • The 
current, whose typical amplitude is twice the value needed to elicit a response for slow 
pacing, activates a propagating excitation wave. Transmembrane voltage is measured with 
a glass microelectrode filled with KC1 conducting fluid. Before collecting any data, the tissue 
is paced at Bq = 1,000 ms for about 30 minutes. We collect data at a sampling rate of 4 
kHz. The data is then processed in custom-written Matlab code. 

We implement the alternate pacing protocol experimentally by repeatedly carrying out 
the following steps: 1) We pace at Bq for two minutes, during which we record the transient 
behavior of the APDs as they reach a steady-state of either a 1:1 or 2:2 rhythm; 2) We 
perturb Bq for 20 seconds at one value of 5 and record the subsequent APDs; 3) We repeat 
step 2 three more times, each time with a new 5 value (for 11 of the 12 trials, the values 
of S sweep from high to low, to include 20 ms, 15 ms, 10 ms and 5 ms; in one trial, we use 
the reverse.); and 4) The perturbations are turned off and pacing at Bq is resumed for 20 
seconds to check that the steady-state value of the system did not drift. The preceding steps 
are repeated for decreasing values of Bq until persistent alternans over several Bqs ensue, or 
the cell fails to respond to every applied stimuli. Twelve of our trials exhibited a transition 
to alternans. Typically, Bq is decreased in steps of 25 ms, which means that the last Bq with 
a 1:1 response is within 25 ms of i?t>if- 

To analyze the results, we determine V vs. 5 at each Bq. In 4 trials from three frogs, 
we find that T shows a decreasing trend as 5 increases, where a typical example is shown in 
Fig. Efa), which agrees with a smooth bifurcation [recall Fig. GJe)]. However, 4 other trials 
from two frogs demonstrate an increasing trend in F as 5 increases, where a typical example 
is shown in Fig. |2fb), which agrees with a border-collision bifurcation [recall Fig. H^f)]. 
In three other trials from three frogs there is no significant variations in T for different <5's 
and therefore cannot be classified into either category. The remaining trial is discussed in 
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FIG. 2: Typical experimental results displaying two different trends in T vs. 5 as revealed by 
alternate pacing for two different frogs. The trend is consistent with (a) a smooth period-doubling 
bifurcation and (b) a border-collision bifurcation. 



detail below. For all experimental data in this paper, the error bars represent the statistical 
and systematic error in the measurement of APD. In the cases where the error bars are not 
evident, the error is smaller than the symbol size. 

Surprisingly, in the remaining experimental trial, we see a smooth behavior when Bq — B^ 
»s 25 ms, but a border collision behavior when Bo — B^f m 50 ms, as shown in Figs. Eft>) 
and (c), respectively. These results are unanticipated because we expected to find evidence 
supporting either one or the other of the two theorized bifurcation types. Figure a) shows 
T vs. B for different values of 5, where it is seen that these curves cross one another. This 
crossing is not indicative of either a smooth or border-collision period-doubling bifurcation, 
and is one of our primary experimental results. 

We find that these seemingly contradictory observations can be explained by a mathe- 
matical model that unfolds a border-collision bifurcation. Consider iterations of the skewed 
tent map (in one dimension) 



(3) 



which is singular at the border x = in state space. This map exhibits a border-collision 
period-doubling bifurcation under the condition: —1 < a+/3 < 1 < a— /? and — 1 < a 2 — (3 2 < 
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FIG. 3: Experimental evidence for both a smooth and border-collision bifurcation in a single 
trial, (a) The trend in T vs. Bq for four different values of 8 (legend). Note that the curves 
cross, indicating the transition from border-collision to smooth dynamics. The bifurcation occurs 
somewhere in the gray area, whose precise location cannot be determined because of our finite 
sampling. The same data is replotted in (b) and (c). The behavior in (b), where (Bq — B^f) < 
25 ms, is consistent with a smooth bifurcation and in (c), where 25 ms < (Bo — B^if) < 50 ms, is 
consistent with a border-collision bifurcation. 
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1. Replacing |x| by a/x 2 + e 2 , where e is a small parameter, results in an unfolding of the 
piecewise-smooth map (jUJ) 

[float fix]x n+ i = /i — a x n — (3 \/x 2 l + e 2 . (4) 

This new map can be viewed as an approximation of the skewed tent map (j3J). In the 
limit e — > 0, map (J3J) reduces to map (j3J). Map Q exhibits a border-collision period- 
doubling bifurcation while, for e / 0, the unfolded map (jlj) exhibits a smooth period- 
doubling bifurcation. However, the dynamics of map (j3J) and map (J3J) exhibit no significant 
differences except when /i is less than or on the order of e. 

We investigate map with a = 2.5, /3 = —2.3, and e = 0.01. The dynamical variable is 
x (similar to APD) and the bifurcation parameter is /i (similar to B Q ). With this nonzero 
e, the bifurcation of map (jlj) occurs earlier than for map (jUJ), specifically at = —0.0091 
rather than for /ibif = for map ©• We apply the alternate pacing protocol to map ()U) by 
perturbing \x with (— l) n 5. Figure Ufa) shows T vs. for different values of 5. These curves 
cross one another [l9j ]. Although the bifurcation is smooth at sufficiently small scales, it 
shows behavior qualitatively similar to a border-collision bifurcation when (jl > 0, in good 
agreement with our experimental findings [see Fig. Of a)]. In Figs. Efb) and (c), V vs. 6 is 
plotted for /i = —0.00773 and [i = 0.01, respectively. For negative values of /U, T decreases as 
5 increases, a trend consistent with a smooth bifurcation [see Fig. El^b)]. On the other hand, 
for positive values of fi, T increases as S increases, a trend consistent with a border-collision 
bifurcation [see Fig. Efc)]. We note that the region over which the dynamics are smooth-like 
[/i < in Fig. El^a)] can be very small and on the order of e in width. Therefore, a smooth 
bifurcation may be masked by a border-collision bifurcation, depending on the proximity to 
the bifurcation point. 

Our results call to question the suggested clinical use of alternate pacing [l^. In a 
typical smooth bifurcation, one expects large Ts far above the bifurcation point. Hence, the 

propensity for alternans could be revealed using alternate pacing for values of Bo far from 

I I 

the bifurcation. It was suggested [15j that such a procedure might be useful in the clinic 
because it would use pacing rates that are slow enough to avoid a potentially life-threatening 
situation. However, we find that T remains small until the pacing rates are decreased to a 
neighborhood of the bifurcation. Consequently, our research suggests that one would not be 
able to probe instabilities clinically through the use of alternate pacing. 
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Our experimental results reveal that the bifurcation, although smooth on a fine scale, may 
exhibit behavior similar to a border-collision bifurcation on a coarse scale. This phenomenon 
is described by a new model called an unfolded border-collision bifurcation. We speculate 
that the origin of this behavior is due to rapid changes in cellular behavior during the course 
of the action potential 20], such as those in the current- voltage or the current-concentration 

n 

relation 211. 
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FIG. 4: Theoretically predicted behavior of the unfolded border-collision bifurcation [map (jl}]. (a) 
r as a function of the bifurcation parameter [i for various values of the perturbation size 5 (legend, 
e = 0.01). The same data is replotted in (b) and (c). T vs. 5 for (b) /j, = —0.0073, a distance of 
0.0018 from the bifurcation point and (c) \i = 0.01, a distance of 0.0191 from the bifurcation point. 
The trend in T vs. 5 is consistent with a smooth bifurcation in case (b) and a border-collision 
bifurcation in case (c). 
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